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We consider the problem of determining the state of a quantum system given one 
or more readings of the expectation value of an observable. The system is assumed 
to be a finite dimensional quantum control system for which we can influence the 
PsJ ■ dynamics by generating all the unitary evolutions in a Lie group. We investigate 

to what extent, by an appropriate sequence of evolutions and measurements, we can 
I obtain information on the initial state of the system. We present a system theoretic 

' viewpoint of this problem in that we study the observability of the system. In this 

, context, we characterize the equivalence classes of indistinguishable states and propose 

algorithms for state identification. 

1 Introduction 

Given a control system 

x = f(t,x,u), (1) 
^ ' 

;_j ■ where u represents the control function, x the state varying on a manifold M, with output 

^ y = y{x), denote by x{t, u, xq) the solution of (1) with control u, initial condition xq, at time 

t. Two states xoi and xo2 are said to be indistinguishable (see e.g. [26]) if, for every control u 
and every time t, we have y{x(t,u,Xoi)) = y{x{t,u,Xo2))- A system is said to be observable 
if no two states in M are indistinguishable. 

In this paper, we are interested in the observability properties of quantum control systems 
whose dynamics are described by the Liouville's equation for the density matrix p (see e.g. 
[3]), 

^J^p=[H{u{t)),p]. (2) 

We shall restrict ourselves to the finite dimensional case where p is an n x n matrix. The 
Hamiltonian H{u{t)) is an n x Hermitian matrix, in general, function of one or more 
control functions u{t). We have assumed here and will assume in the rest of the paper 
that we are dealing with closed (noninteracting with the environment if not through the 
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control functions and during the measurement process) quantum system. We assume that 
we perform a measurement of the mean value of an observable, represented by a Hermitian 
matrix S. In this case the output y takes the form 

y^TriSp). (3) 

Since Tr{p) = 1, it will be convenient to replace p with the traceless matrix p — ^Inxn and 
S with the traceless matrix S — '^-^-^Inxn- This has the effect of 'shifting' the value of the 
output by a constant Tr{S) value which does not play any role in the indistinguishability 
considerations that will follow. Therefore we will set in the following Tr(p) = and Tr{S) = 
0. The solution of (2) varies as 

p{t)=X{t)p{0)X*{t), (4) 
with X solution of the Schrodinger equation, 

X{t) ^ -iH{u{t))X, X(0) = 4xn. (5) 

Prom known results in the theory of quantum control (see e.g. [16], [19], [23], and see [9], [18] 
for the non-bilincar case), X can be driven to every value in the Lie group corresponding 
to the Lie algebra C generated by spanu^y{iH {u)} where U denotes the set of possible values 
for the control. With initial condition /Ci(0), Hermitian and with trace zero, the density matrix 
p can attain all the values in the orbit 

O := {Xp{Q)X*\X e e^}. (6) 

A study of the observability of control systems involves two main things. First, one would 
like to collect, in equivalence classes, initial states that cannot be distinguished by varying 
the control and measuring the output (see the next two sections for definitions in our case) . 
Second, one would like to have methods to infer the equivalence class of the initial state 
from appropriate sequences of measurements and evolutions. We consider these problems 
for quantum control systems in this paper. 

The question of determination of the state of a quantum system from measurements is 
at the heart of quantum mechanics and it was already discussed by Pauli in [22]. Several 
contributions have appeared in recent years and a discussion of the problem in general terms 
can be found in [6], where, like in the present paper, the problem of determination of the 
initial state (as opposed to the current state) was described. We present in this paper a 
treatment of this topic from a system theoretic view point. In this context, our study is 
closely related to other studies of the observability of nonlinear systems [15], [17], [21] (see 
also [8] for systems varying on Lie groups). However we consider here a specific model 
for which we can obtain more complete results. Moreover a new element appears in the 
treatment of quantum systems, that is the transformation of the state as a result of each 
measurement. This can take different forms according to the type of measurement considered 
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(see e.g. [4], [7], [10]). We shall mainly consider the case of Von Neumann measurement 
[24] , [27] and discuss extensions to other cases. 

The paper is organized as follows. In Section 2 we define and describe the set of states 
that cannot be distinguished in one measurement. Then, we generalize in Section 3 to states 
that cannot be distinguished in multiple measurements. The determination of the state from 
one or more measurements is discussed in Section 4. Conclusions are given in Section 5. 



2 Indistinguishability and observability with a single 
measurement 

In the following, S is the traceless Hermitian matrix representing the observable and p{t, u, po) 
is the solution of (2) at time t, with initial condition equal to po, and control u. 

Definition 1. Two states pi and p2 are indistinguishable in one step if, for every control 
function(s) u and every t, we have 

Tr{Spit, u, pi)) = TriSpit, u, P2)). (7) 



The definition asserts that two states pi and p2 are indistinguishable if there is no ad- 
missible experiment involving only one measurement which would give different results with 
initial states pi and p2- It is clear that indistinguishability in one step is an equivalence 
relation. The set of possible values for the density matrix will be denoted by TZ. It is a 
convex subset of the vector space of n x n Hermitian matrices (with zero trace), isu{n) and, 
in general, the vector space spanned by the elements of TZ is the same as i{su{n)). The 
elements of TZ are parametrized by — 1 = dim su{n) parameters ^ 

Definition 2. The system is observable in one step if pi and p2 & TZ are indistinguishable 
in one step only when pi — P2- 

Instrumental in the characterization of classes of indistinguishable states is the vector 
space oi n X n skew- Hermitian matrices, 

V - ®T=o(^d1iS. (8) 

Here ad)4S is the space obtained by taking k Lie brackets of iS with elements in the Lie 
algebra C We shall call V, Observability Space. If is a set of generators of the Lie 

algebra C, it follows from an application of the Jacobi identity (see Appendix A) that the 
observability space V is spanned by the matrices^ 

ad)} ad'}? ■ ■ ■ ad'}? iS, (9) 

^In the presence of special symmetries, a parametrization with fewer parameters can be given (see [2] for 
an example). 

^ad^T := [R, [R, ...[R,T]]] where the Lie bracket is taken k times. 
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with ki,...,kr > 0, and {ji, ■••,>} G {l,...,m}. V is the smallest subspace of su{n) stable 
under C ^ and containing iS. V might not be a Lie Algebra, however it is always a subspace 
of the Lie Algebra generated hj iS, Bi, and therefore a subspace of su{n). Therefore 
its dimension is bounded by dim. su{n) = — 1. Notice that V can be calculated with an 
algorithm that, at each step, calculates the matrices of 'depth' d + \ from the matrices of 
depth (i, where the depth is the number of Lie brackets performed namely ki -\- k2 -\- ■ ■ ■ -\- k^- 
in (9). The algorithm starts with the matrix iS, which has depth 0, and ends when the 
dimension reaches — 1 or there is no increment in the dimension. By finite dimensionality, 
there is always a finite k such that 

V = ^l^^adlliS. (10) 

We have the following result that relates the partition of the state space into classes of 
indistinguishable states with the properties of the observability space V. 

Theorem 1. The following three conditions are equivalent 

1. The states pi and p2 are indistinguishable in one step. 

2. For every X & e^, 

Tr{X*SXpi) = Tr{X*SXp2). (11) 

3. For every F e V, 

Tr{Fp,)^Tr{Fp2). (12) 



Proof. The equivalence between conditions 1 and 2 simply follows from the fact that the set 
of values obtainable for p starting from p{0) = pi^2 is described in (6), and from elementary 
properties of the trace. 

Now assume (11) holds and choose k matrices Ri,...,Rk (not necessarily all different) in 
jC. Then, for every /c— pie of real numbers ti, ...,tk we have 

Tr(e~-^i*i • • • e-^'^^HSe^''*'' ■ ■ ■ e^'^'pi) = Tr^e"^^*^ ■ ■ ■ e'^'^^HSe^'^*'' ■ ■ ■ e^^^Vz)- (13) 

Calculating the derivative, , of both sides we obtain 

Tr{adR^adR^ ■ ■ ■ adn^Spi) = TriadR^adn^ ■ ■ ■ adRjSp2), (14) 

which proves Condition 3, since k and Rj, j — l,...k, are not specified. To prove that 
Condition 3 implies Condition 1, let Fi, be a basis of V with Fi = iS. Then we have, 
using (2), 

jTr{Fjp{t,u,p,,2)) = E aj,k{t)Tr{Fkp{t,u,p,,2)). (15) 
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for some (time varying) coefficients aj,kif) depending on the control u. Therefore we have 
that Tr{Fjp{t,u, p-i)) and Tr{Fjp{t,u, p2)), satisfy the same (hnear) system of differential 



equations and since the initial conditions are the same, then 

Tr{F,p{t, u, pi)) = Tr{Fjp{t, u, pa)), j = 1, s. (16) 
In particular, we have, 

Tr[Sp{tAi,p^)) = Tr{Sp{t,u,p2)). (17) 
Therefore the two states are indistinguishable. □ 



The inner product < •,• > in su{n) is defined as < A, 5 >= Tr{AB*). Theorem 1 
states that two matrices in are indistinguishable if and only if they differ by an element 
in V"*". Therefore we can state the following criterion of observability in one step which is a 
consequence of Theorem 1. 

Theorem 2. System (2) is observable in one step if and only if one of the following equivalent 
conditions are verified 

1. 

spanxee^X*iSX — su{n), (18) 

2. 

V^su{n). (19) 

Remark: The notion of observability is closely related to the notion of informational com- 
pleteness of observables as treated for example in [14]. A set of observablcs B is called 
informationally complete if Tr{Bpi) = Tr{Bp2) for every B E B implies pi = p2- From 
condition 2 of Theorem 1 and the definition of observability, we can say that a system is 
observable if and only if the set of operators {X*SX\X e e-^} is informationally complete. 

2.1 Relation between controllability and observability in one step 

\i C = su{n), namely the system is operator controllable [1], and iS 7^ 0, then it is also 
observable. In fact, in this case, we have 

spanxee^X*iSX — spanxesu{n)X* X — su{n). (20) 

To verify this we can more easily verify condition (19). Since V is a nonzero ideal of su{n) 
and su{n) is a simple Lie algebra, V must be equal to su{n). Therefore we have. 

Corollary 3 Controllability along with S ^ implies observability in one step. 

The converse of Corollary 3 is not true not only because we may have the equality 

{X*SX\X e e^} = {X*SX\X e SU{n)}, (21) 
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even though £ ^ su{n) [1], [25] but also because we may have (18) (19) even though (21) is 
not verified. A simple example of this can be found already in the n — 2 case by taking 

_^.), £^span{(^\ J)}. (22) 

2.2 First order conditions for observability in one step. 

The case of the equality of the orbits in (21) is particularly favorable because we can give 
a different condition of observability which avoids the calculation of repeated Lie brackets 
for V and involves only the calculation of Lie brackets of depth 1. We have the following 
Proposition. 

Proposition 4 Assume S 0. The system is observable in one step if one of the following 
two equivalent conditions is verified 

1. 

{X*SX\X e e^} = {X*SX\X e 5C/(n)}, (23) 

2. 

[C,iS] = [su{n),iS]. (24) 

Prom a practical point of view condition (24) may be easier to verify since it involves 
calculation of first order Lie brackets only. The condition tells us that, by calculating first 
order Lie brackets, wc can infer the properties of V which is defined through higher order 
Lie brackets. If condition (24) is not verified we may still have observability. 

Proof. The equivalence between the conditions (20) and (23) was proven in [1] although in 
a different context"^, therefore we shall not repeat the proof here. Clearly (23) implies (18) 
with (20) and (19) and therefore observability. □ 

Condition (24) can be verified by comparing the dimensions of the two vector spaces. The 
dimension of [iS, su{n)\ can be expressed in terms of the multiplicity of the eigenvalues of iS 
(recall that iS is not zero and it has zero trace so it has at least two distinct eigenvalues). 
We have 

dim[iS', su{n)] = 2 ^ ^j^fe, (25) 

3<k 

where rij {rik) is the multiplicity of the j— th (A;— th) eigenvalue. 

If iS is known to be in a proper subspace T of su{n) stable under C (e.g. C itself or C^) 
then we cannot have observability because V T ^ su{n). 

^There, S was the density matrix and we wanted to give practical conditions to verify that the set of 
possible density matrices that can be obtained by varying X in e"^ is the same as the largest possible one 
namely the one obtained by varying X € SU{n). This condition was called Density Matrix Controllability. 
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Example: Two spin | particles are interacting through Ising interaction and are driven by 
an electro-magnetic field in the x direction [11], [12]. The magnetic field couples with one of 
the spins only and we can detect the magnetization in the z direction. Denote by (7x,y,z the 
x,y,z Pauli matrices (see e.g. [24]) 

^^^=K? J)' ^^^=K^ 7)' ''''■=lCo -l)' ^^^^ 

by 1 the 2x2 identity matrix and by u — u{t) the x component of magnetic field. After 
appropriately scaling the parameters involved, the Hamiltonian H has the form 

H = a^i^az + u{t)axiS)l, (27) 

and the output matrix S is given by = cr^ 1 + 1 ® a^. The dynamical Lie algebra C is 
spanned by icx^ Cz, H ® cTx and ia^ (8> Cy. We have from formula (25) dimfiS*, su{n)] = 8 
while dim[zS', C] = 2. Therefore the sufficient criterion of observability of Proposition 4 fails. 
Moreover since iS is in and JC^ is stable under £, the system is not observable. 

2.3 Decomposition of the state space 

It is natural to decompose the state p as p{t) — pi{t) + p2{t), with pi{t) e V and p2{t) e V"*", 
for every t. Then we have 

Pi = -^[i7(w),pi], (28) 
P2^-i[H{u),p2l (29) 

and 

TT{Sp{t))=TT{Sp,{t)), (30) 

for every t. Therefore if we are interested in the effect on the output S we can parametrize 
only the component of /9 in V. 

3 Indistinguishability and observability with multiple 
measurements 

We now generalize the above characterization of states that are indistinguishable after one 
measurement to states that are indistinguishable after k measurements, for general k. In fact 
it may happen that, even if two states give the same output function at the first measurement, 
for every control and at every time, they give different values at the second measurement. 
This is a consequence of the fact that the first measurement modifies the state. Modern 
quantum measurement theory (see e.g. [4], [7], [10]) has studied ways to model the change 
in the state due to measurement as well as ways to integrate the measurement process in 
the framework of quantum dynamics. We shall remark at the end of this section on possible 
extensions to other types of measurements but will consider the simplest case where the 
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quantum measurement postulate [24] [27] holds. This is also called Von Neumann (or Von 
Neumann-Liiders) measurement. More precisely, rewrite the observable matrix S as 

n n 

S=Y. \^3^] H \^3^ (31) 
i=i i=i 

where aj are the orthonormal eigenvectors of 5", ITj, j = l,...,n, are the associated pro- 
jection matrices defined by Uj := aja* and are the associated eigenvalues. Define the 
automorphism V in the space of (skew)Hermitian matrices 

n 

V{F)^J2^jFUj, (32) 

which rctm^ns the diagonal part of F, if wc arc working in a basis where 5* is diagonal. If the 
state at the time of the measurement is p{t, u, po), according to the measurement postulate, 
the state after the measurement is V{p{t,u, po)). Assume the experiment consists of an 
evolution for time ti with control Ui, followed by a measurement, followed by an evolution 
for time t2 with control U2, followed by a measurement, and so on, up to an evolution for 
time tk with control Uk- The k—th measurement, at time ti + 12 + ■ ■ -tk, gives the result 

yk{ti, tk, Ml, Uk, po) ■■= Tr{Sp{tk, Uk, V{p{tk-i, Uk-i, V{- ■ ■V{p{ti,ui, po))))))). (33) 

We can extend Definitions 1 and 2 as follows 

Definition 3. Two states pi and p2 arc indistinguishable in k steps if for every sequence of 
control function(s), Ui,U2, ■■■,Uk, defined in intervals [0, ti), [0, ^2) v)[0, ^fe]; we have 

yk{ti, ...,tk,ui, ...,Uk,pi) = Vkih, ...,tk,ui, ...,Uk,p2), (34) 

Definition 4. A system is observable in k steps if no two states are indistinguishable in k 
steps. 

Definition 5 Two states are indistinguishable if they are indistinguishable in k steps for 
every A; > 1. A system is said to be observable if no two states are indistinguishable. 

It is convenient to rewrite the output at the A;— th measurement in terms of the values 
of the evolution operator X in (5) at the endpoints of the intervals [0, ti),...,[0, t^]. We call 
these values of X, Xi,...,Xk. Using (4), we have 

Vk := VkiX,, ...,Xk,po) = Tr{SXkV{Xk-iV{- ■ -ViX^poXl) • ■■)Xl_,)Xk). (35) 

Therefore an alternative definition of indistinguishability in k steps can be given, that is pi 
and P2 are indistinguishable if for every set of values Xi,...,Xk in e^, yk{Xi, Xk, pi) — 
yk(Xi, ...,Xk,p2)- 

We can give conditions of indistinguishabihty and observability as in Theorems 1 and 2, 
by introducing Generalized Observability Spaces. More specifically, define the Observability 
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Space of order 0, Vo := span{iS}, and the Observability Space of order 1, Vi := V in (8). 
The Observability Space of Order k, Vk, is defined recursively by 

Vk:=®f=oadiV{Vk-i). (36) 

It is the largest subspace of su{n) containing V{Vk-i) and stable under C It also follows 
from a proof analogous to the one in Appendix A that, if is a set of generators of 

the Lie algebra £, is spanned by the matrices 

ad%y%^---ad%^F, (37) 

with F e V{Vk^i), ki, kj. > 0, and {j2, Jr} G {1, ■■■,rn}. Notice also that it follows by 
induction, since Vq QVi, that 

Vfc-i C Vk, (38) 

for every A; > 1. 

We have the following generalization of Theorem 1. 
Theorem 5. The following three conditions are equivalent 

1. The states pi and p2 are indistinguishable in k steps. 

2. For every k—ple Xi, ...,Xk with values in e^, 

Tr{x;v{x*v{- ■ ■r{xi_,r{xisXk)Xk-i) ■ ■■)X2)x,p^) = (39) 

Tr{XlV{X*V{- ■ ■V{Xl_{P{XlSXk)Xk-i) ■ ■■)X2)X,p2). 

3. For every F eVk, 

Tr{Fp,)^Tr{Fp2). (40) 



It follows from (40) and (38) that if two states are indistinguishable in k steps they are 
indistinguishable in r steps for every r < k. In other terms if we can distinguish two states 
in r steps we can distinguish them in A; > r steps as well. 

Proof: If pi and p2 are indistinguishable, then, for all the Xi,...,Xk in e^, we have 
VkiXi, Xk, pi) = yfe(Xi, Xfc, P2) in (35). Now notice that, for a general po, 

TT{sXkr{Xk-{P{- ■ ■vix.poxi) . ■■)xi_,)x;) = (4i) 
Tr(x;sXkV(Xk-iV(- ■ ■vix.poxi) ■ ■■)x;_,)) = 
Tr{v{x;sXk)Xk-iV{- ■ -vix^poxi) ■ ■■)x;_,) = 

Tr{X;_,V{X;SXk)Xk^^V(- ■ -ViX.p^Xl) ■ ■■)) = 



Tr{xiv{x;v{- ■ ■v{x;_,v{x;sXk)Xk-i) ■ ••)^2)^ipo). 
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Using this for po = Pi and po = P2 along with (35) we see that indistinguishability of pi 
and P2 in k steps implies equation (39). The proof that Condition 2. imphes Condition 3. 
is exactly analogous to the corresponding proof in Theorem 1. The proof that Condition 3. 
implies indistinguishability also is a generalization of the corresponding proof in Theorem 1, 
with some more elements that we now illustrate. Consider a basis Fj, j = 1, s, of Vfc and 
derive a differential equation for Tr{Fjp{t, Mi, pi,2))- The differential equations corresponding 
to pi and P2 are the same with the same initial conditions, because of the assumption (40). 
Therefore, in particular, at time ti, we have 

Tr{Fjp{t,,u,,p,)) = Tr{Fjp{t,,u,,p2)), (42) 

for every Fj e Vk and therefore for every Fj e V{Vk-i)- If Fj, j = 1, s, is a basis of Vk-i, 
then we have 

Tr{V{Fj)p{ti,m,pi)) = Tr{V{Fj)p{ti,m,p2)), (43) 

that is 

Tr{F,V{p{t,,uuPi)))=Tr{FjV{piti,ui,p2))). (44) 

Now, derive a differential equation for the variables Tr{Fjp), with Fj a basis of Vk-i, on the 
second interval of length t2 and with control U2. The function corresponding to pi satisfy the 
same differential equation as the function corresponding to p2 and since the initial conditions 
are the same, from (44), we obtain that for every Fj in Vk-i 

Tr{Fjp{t2, U2, P(p(ti, Ml, Pi)))) = Tr{Fjp{t2, U2, V{p{h,u,, ps)))). (45) 

This is, in particular, true for elements of V{Vk-2)- Proceeding this way, after k steps, we 
obtain the equalities of outputs i/k in (33) for po = pi and po = P2, for every /c— tuple ti, ...,tk 
and controls Ui, ...,Uk, and therefore indistinguishability. □ 



An example of Vi 7^ V2 is given by 

/I 0\ 

S :^ 0-3 
VO 2/ 



jC :— span{ 



/ i 2\ 

-i 

\-2 oy 



}• 



(46) 



We also have the following Theorem concerning observability. 



Theorem 6. System (2) is observable in k steps if and only if one of the following equivalent 
conditions is verified 



1. 



2. 



spanx,,x,,...,x,eecX;V{X;r{- ■ ■V{Xl_,V{XliSXk)Xk-,) ■ ■■)X2)X, = su{n), (47) 

Vk^su{n). (48) 

A system is observable if and only if there exists a k such that one of the equivalent 
conditions (^7), (48) is verified. 
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To check observability we only need to verify (48) for a finite number of k^s until we find 
a k such that Vk-i — Vk or Vk — su{n). 

It is obvious that since controllability {C — su{n)) implies observability in one step it 
also implies observabihty in k steps for every k. The natural extension of the condition (23) 
of Proposition 4 would be 

{x*v{x;v{- ■ ■v{x;_,v{x;sXk)Xk-i) ■ ■■)X2)x,\x,,X2, ...,Xk e e^} = (49) 
{xiv{x;v{- ■ ■v{x;_,vixisXk)Xk-i) ■ ■■)X2)x,\x,,X2,...,x, e su{n)}. 

However we cannot give a Lie Algebraic condition for (49) (which would be an extension of 
(24) for this case). Notice that (24) is essentially the equality of the tangent spaces at S of 
the two manifolds in (23). The main difficulty is that the two sets in (49) are not guaranteed 
to be manifolds. For example, if we consider 




and := SO{2), and k = 2, then we have 

{X*P(X*5X2)Xi|Xi,X2ee^} = {(^ ^^^\a,be R,V^T¥ <1}, (51) 

which is a manifold with boundary. 

Like for the case of indistinguishability in 1 step, we can write 

p{t)^p,{t) + p2{t), (52) 

with pi{t) e Vk and p2{t) e Vj^, which satisfy the equations (28), (29), (30). Therefore if 
we are interested in the effect on the output S we can parametrize only the component of 
p in Vk- In particular, if Vk is the largest of the observability spaces we can neglect the 
component p2{t) of the state since it will not have any effect on any measurement. 

Remark: The above treatment, which has been presented for Von Neumann measurements, 
can be extended to more general types of measurements (see e.g. [4], [7], [5], [10], [13]). 
We have used the fact that, according to the measurement postulate, the state changes as 
p — > ■P(p). For a more general measurement, with a countable set of possible outcomes M., 
the state will change according to 

P^Hp)--= E ^m{p). (53) 

The super-operators are called operations and according to Kraus' representation theorem 
[20] , under suitable assumptions, can be expressed as 

k 

11 



for a countable set of operators ^mk- Our treatment will go through by replacing V{p) with 
^(p). In particular, we can define a dual super-operator J^* acting on observables as 

J'*{s):= J2 <^>i{s), <^*m{s):=J2ni,snr,,. (55) 

mEM. k 

This has the property Tr{J^*{S)p) = Tr{SJ^{p)) and we can use this to extend the calcula- 
tions in Theorem 5. Moreover the definition of Vfc in (36) has to be replaced by 

Vfc:=©°^o«4-^*(Vfc-i). (56) 



4 Initial State Determination 

We now investigate how much information on the initial state we can extract from an ex- 
periment which alternates prescribed evolutions with measurements. We deal with a single 
experiment and with a single quantum system rather than with many copies of the same 
system, as it is done some times in this context. We shall assume, for simplicity, that the 
system is controllable namely C — su{n). Moreover, we can assume, without loss of general- 
ity, that the output matrix S is diagonal. We shall use the following formula (see (35), (39)) 
for the output at the k-ih. measurement 

Vk = Tr{X,poX*V{X;V{- ■ ■V{Xl_,r{XlSXk)Xk^i) ■ ■■)X2)), (57) 

for the unknown initial state po- Now, since every matrix of the type V{-) is diagonal, it 
follows from (57) that it is only possible to obtain information on the diagonal elements of 
XiPqX^ and therefore on at most n — 1 independent parameters of the unknown matrix po- 
It is in fact possible to obtain all the n — 1 independent diagonal elements of the matrix 
XiPqXi :— Po- At the first measurement we obtain 

yi = Tr{poS). (58) 

Then we choose X2 as a permutation matrix so that S2 '■= X2SX2 is still diagonal but the 
diagonal elements are a permutation of the diagonal elements of S. We also have V{S2) = S2 
so that, at the second measurement, we obtain 

y2 = Tr{poS2). (59) 

Then we choose the evolution X3 with X3 := X3X2 and X3 performing another permutation 
of the diagonal elements of S. X^SX^ :— S^. Therefore, the third measurement gives 

ys^TripoSs). (60) 

Continuing this way, we can obtain n\ equations for the diagonal elements of po, Xi, ...,Xn, 
i.e. 

n 

J2 ajkXk = yj, J = 1, n!, (61) 

k=l 
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where the elements ajk are appropriate permutations of the diagonal elements of S. To this 
we have to add the equation^ 

7^KPo) = Ea;fc = 0. (62) 

fc=i 

If S is not a scalar matrix, it is always possible to choose n — 1 permutations and therefore 
n — 1 equations in (61) that together with (62) have a unique solution. In fact the n! + 1 x n 
matrix obtained by placing in the first n\ rows all the permutations of the diagonal elements 
of S and in the last row 1, 1, 1 has always rank n. The rank of this matrix is the same 
as the rank of a matrix obtained by adding to every row the last row (1, 1, 1) multiplied 
by an arbitrary constant. Therefore we can assume that the elements of the matrix are 
nonnegative and apply a Lemma in Appendix B. 

As seen above, in the Von Neumann case, the number of independent parameters that 
can be inferred by a sequence of evolutions and measurements is bounded by the dimension 
of the system. This suggests to consider different types of measurements to obtain complete 
information on the initial state of the system. One possible scheme is as follows. Consider 
a system Ei of dimension n, with unknown state pi and couple it with a (large) system E2, 
of dimension m, whose state is known to be p2- The density matrix of the coupled system p 
at time is 

p(0) = pi®P2. (63) 

This matrix has dimension nm and only — 1 parameters are not known. Now, if we let p 
evolve, after time the matrix p{t) cannot in general be written as a tensor product, since 
the two systems are now entangled [28]. If we perform repeated Von Neumann measurement 
on the coupled system, we are able to obtain information on nm — 1 independent parameters 
of p. Since p contains — 1 unknown parameters only, we may be able to obtain information 
on all of them if m > n. We give now a simple numerical example of this scheme. 

The unknown state of a spin | particle is represented by the density matrix (without 
shift of the trace) 

m I 



with m real. Two spin | particles with known state 

i 0\ /I 



are coupled with it. Therefore the unknown state 

Po Pi ® P2 (66) 

has only three unknown parameters. We can observe the magnetization in the z direction 
of the system of two spins. The associated matrix is given by 

S^a^<S>l<S>l + l(S>(T^(S>l + l<^l<S><T^ (67) 

^Recall that, without loss of generality, we are considering density matrices with trace equal to zero rather 
than one 
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(see (26)) which is diagonal. Prom formula (57) and the previous discussion we can obtain 
the diagonal elements of the matrix XiPqX^. Consider the vectors 



If the first three columns of are chosen as 



xs := (g) 62 <8) Ci) 
then we obtain for the diagonal elements 



9 

^2Po^2 = ^(l + 2i?e(/)) 

x*2poX3 = l{l + 2Im{l)). 



9 

Prom this we can extract the values of m and I 



Xi :— Ci (8) ei Ci (69) 
X2 := vi<Siei<Si e-i 



x{pqX\ = -m (70) 



5 Discussion and Conclusion 

In this paper, we have presented a treatment of the observabihty properties of quantum 
systems compatible with quantum measurement theory. We have focused on Von Neumann 
measurements but indicated extensions to more general types of measurements. We have 
given a characterization of states that cannot be distinguished in one or more measurements 
and conditions for observabihty. Contrary to most studies on observabihty of nonhnear 
systems (see e.g. [26]) conditions for observabihty and indistinguishabihty are global in this 
case, however observability does not always imply that it is possible to infer from appropriate 
evolutions and measurements all the parameters of the initial state. In fact, for Von Neumann 
measurements, there is a natural limit to the number of parameters of the state that can be 
derived. This does not improve if we consider measurements of different (and not necessarily 
commuting) observables. In this case, the resuh of the k—th measurement has still the form 
(57) although now each projection V corresponds to a possibly different observable and S 
corresponds to the observable measured last. In this case the first V on the left is always the 
projection corresponding to the first measurement and therefore, once again, only at most 
n — 1 independent parameters of XiPqX^ can be obtained. We have seen, in the previous 
section, that complete information on the initial state may be obtained by coupling the 
system with an auxiliary system whose state is known. 
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Appendix A: Evaluation of V using a set of generators 
of jC. 

Let Bi, Bm a set of generators of C. Denote by V the space spanned by the matrices in 
(9). It is obvious that 

V C e~ oasis'. (71) 

To show 

®^=^adliS C V, (72) 
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we first show that 

[V,>C]CV. (73) 

It is enough to show for elements F in a basis of L given by Bi, -B^, and hnearly inde- 
pendent (repeated) Lie brackets, [F, V] C V. We proceed by induction on the depth of F . 
If F is of depth 0, namely F is one of the matrices -Bi, ...,i?m then (73) follows from the 
definition of V. Now, let us assume (73) true for matrices F of depth < d and let us show 
it for matrices F of depth d + 1. In particular, write F q& F :— [Z, T], where Z is of depth 
d and T is of depth zero. If F is a matrix in V, from the Jacobi identity, we obtain 

[V, [Z, T]] = -[Z, [T, V]] - [T, [V, Z]], (74) 

since both terms on the right hand side are in V, from the inductive assumption, we have 
that the term on the left hand side is also in V, therefore we have proved (73). Now, from 
(73) we have 

adciS ■.= [iS,C]CV, (75) 

and from this 

adliS := [[iS, £], £] C [V, £] C V, (76) 
where we have used (73). Proceeding this way, we see that for every k > 

ad'liS C V, (77) 

which proves (72). 



Appendix B 

Lemma Let he n non-negative nvmbers not all equal. Consider the matrix 

A{xi, . . . ,Xn) G j^"'^^" whose rows are the the permutations of xi, ...,Xn- Then the matrix 
A = A{xi, . . . , Xn) has rank n. 
Proof 

Let 2 < r < n be the number of different values assumed by the {a;j}i=i ... „. Denote by 
< di < ■ ■ ■ < dr these values, and let li, for z = 1, . . . , r be the cardinality of { j | Xj — d^ }. 
Thus h = n. We will prove our statement on induction on r > 2. 

case r = 2 

We prove this part by induction on n > 2. If n = 2, then the statement is easily proved 
by computing the determinant of A. Let n > 2. Since all the columns of A sum up to the 
same value, which is strictly positive, setting 




wc liavc that 

rank A = rank A' = rank 
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where we have set 

Choose a value xj e {xi, . . . ,x„}, such that xj — di, assuming that h > 2 (otherwise 
choose it so that xj = ^2)- Assume that we have rearranged the rows of A' in such a way 
that the first element of the second to the (n — 1)! + 1-th row is x^. Notice that this can be 
done since the rank remains unchanged. Then for z = 2, . . . , (n — 1)! + 1 we subtract from 
the i-row of A' the first row multiplied by xj. Notice that ii di — we leave the matrix 
unchanged. After this operation, the matrix A' has the following form: 



/ 1 1,---,1 \ 



A'^ 



V 



B 



(80) 



Notice that A is an (n — 1)! x (n — l)-matrix with the same structure of A and values 
Vj ~ -^j particular the y/s are either or ^2 — <ii. Thus, by inductive assumption 

we have that rank A — n — 1, which, in turn, implies rank A' — n as desired. Had we 
chosen xj = d2, we would have had all the values yj < with the two possible values and 
di — d2, then we would have changed the sign of A (which does not affect the rank) and 
applied the inductive assumption. 

case r > 2 We assume that the result is true for r — 1. The idea of the proof is similar to 
the r — 2 case. Assume again that we have chosen xj e {xi, . . . , Xn}, such that xj — di and 
we have performed to the matrix A' (defined in (78)) the same operation as in the previous 
case to put A' in the form (80). 

Now, if we prove that rank A — n — 1 then we get rank A' = n. As before, A is an 
(n — 1)! X (n — l)-matrix which is the same structure as A with values yj = Xj — xj > for 
j ^ i. If Zi = 1, then we are done by the inductive assumption since the numbers yj assume 
r — 1 different nonnegative values. If /i > 1, then we perform the same procedure as before 
starting with A instead of A. Notice that A has n — 1 different numbers, and is such that 
di = and the cardinality of the {yj = 0} is h — 1. Thus we need to repeat this procedure 
li times and then we can conclude by induction. 

□ 
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